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ABSTRACT: A methodology is shown for predicting the time-dependent reliability 
of ceramic components against catastrophic rupture when subjected to transient 
thermomechanical loads (including cyclic loads). The methodology takes into 
account the changes in material response that can occur with temperature or time 
(i.e., changing fatigue and Weibull parameters with temperature or time). This 
capability has been added to the NASA CARES/Z//e (Ceramic Analysis and 
Reliability Evaluation of Structures/ Z.//T) code. The code has been modified to have 
the ability to interface with commercially available finite element analysis (FEA) 
codes executed for transient load histories. Examples are provided to demonstrate the 
features of the methodology as implemented in the CARES/Z//e program. 

Nomenclature 

A slow crack growth coefficient (Ai - power law, A2 - Walker law) 

A e effective area 

a crack length 

B slow crack growth material parameter 

C Shetty shear sensitivity constant 

CARES Ceramics Analysis and Reliability Evaluation of Structures 

f c cyclic frequency 

FEA Finite Element Analysis 

Kj Stress intensity factor (i = I, II, III for modes I, II, or III, respectively) 
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Kic critical mode I stress intensity factor (fracture toughness) 

K leq equivalent mode I stress intensity factor from applied effective stress 
Kjcq cy C iic max maximum mode I equivalent stress-intensity factor over a cycle 

k B normalized Batdorf crack density coefficient 

m Weibull modulus (scatter parameter) 

N slow crack growth exponent 

Pf a attenuated probability of failure from proof test and service load 

Pf probability of failure (Pf=l- P s ) 

P s (tf) probability of survival at time t = tf 

Psa attenuated probability of survival from proof test and service load 

Psi probability of survival from proof test and service load 

Ps p probability of survival from proof test 

Q Walker law R-ratio sensitivity exponent 

R R-ratio; ratio of minimum stress divided by maximum cyclic stress 

S surface area based property (indicates surface flaw analysis) 

SCG Slow Crack Growth 

T temperature 

t time 

t p proof test time 

t q total time in proof testing and service 

V volume, or volume based property (indicates volume flaw analysis) 

V e effective volume 

Y crack geometry correction factor 

x,y,z represents a location in the body of the structure 

Z number of load blocks or number of cycles 

A increment or difference 

Atj duration of time step j 

<jf represents the peak stress in a component 

Gieq equivalent mode I far-field stress 

a leq.j equivalent mode I far-field stress during time step j 

^ieq,k,Tmax maximum transfonned stress during k time steps 

o n applied far-field stress nonnal to a crack face 

o 0 Weibull scale parameter 

oob Weibull scale parameter incorporating the effect of k B 

oe characteristic strength (value of Of where 63.21% of specimens fail) 

t applied far-field shear stress on a crack face 

'F represents a location (x,y,z) and crack orientation (a,P) 

dO sin a da dp 

Introduction 

Ceramic materials display a stochastic strength response because of the combination 
of brittleness and the random nature of the size, orientation, and distribution of inherent 
microscopic flaws. In addition, the ability of ceramic structures to sustain loads degrades 
over time due to a variety of effects such as slow crack growth (SCG), cyclic fatigue, 
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creep, and oxidation. Life prediction codes such as CARES/Zz/e [/], CERAMIC/ERICA 
[2], and STAU [3] are available and have been demonstrated to be successful in 
predicting the probability of ceramic components failing from spontaneous catastrophic 
rupture when these components are subjected to multiaxial loading and SCG conditions. 
The objective of this paper is to describe and demonstrate the methodology employed in 
the CARES/ Life code that enables the calculation of component failure probability when 
loading and temperature varies over time. This capability is referred to as transient 
reliability analysis and can be used to predict component reliability for situations such as 
thennal shock, startup and shutdown conditions in heat engines, and cyclic loading. This 
computational methodology is an extension to the work of Paluszny and Nichols [4], 
Stanley and Chau [5] and Ziegler [6]. 

For the CARES/ Life program the transient reliability analysis methodology is 
developed with the following features: 

a) Fast-fracture transient analysis (reliability analysis without SCG). 

b) SCG transient modeling with SCG and Weibull parameters (including the shape 
parameter, m) allowed to vary with time and temperature. 

c) Computationally efficient algorithm to compute reliability for cyclic loading. 

d) Cyclic fatigue modeling using a combined SCG and Walker fatigue law. 

e) Transient proof test capability. 

Three example problems are shown; 1) A disk in thennal shock to illustrate fast- 
fracture transient analysis, 2) a diesel engine exhaust valve to contrast predictions for 
cyclic loading and proof testing, and 3) the effect of changing Weibull and fatigue 
parameters over time on the reliability of alumina flexure bars in static fatigue. 

Methodology 

The detailed development of the transient reliability methodology for SCG for the 
CARES/Zz/e program has been described in [7,5] and will not be repeated here. Instead a 
highlight of these equations will be presented as well as a proposed relation for cyclic 
fatigue. 

The CARES/Zz'/e software describes the probabilistic nature of material strength using the 
Weibull cumulative distribution function \9\. For uniaxially stressed components the 2- 
parameter Weibull distribution for volume residing flaws describes the component fast- 
fracture failure probability, Pfv, as 


Pfv = 1 - ex P 


1 


'ov 


r 

G 
J v 


(x,y,z) mv dV 


( 1 ) 


where V is the volume, terms that are a function of volume have V in the subscript, 
a(x,y,z) is the uniaxial stress at a point location in the body, and m and Go are the shape 
and scale parameters of the Weibull distribution, respectively. The shape parameter is a 
(unitless) measure of the dispersion of strength while the scale parameter is the characteristic 
strength of a unit volume of material in uniaxial tension and has units of stress-volume 1 m . An 
analogous equation based on surface area can be shown for flaws that exclusively reside on 
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the component surface. However, for the sake of brevity this paper is restricted to showing 
volume flaw based relations only. 

Determination of the Weibull parameters comes from rupture experiments of specimens 
in simple tension or flexure. Regression techniques such as least squares and maximum 
likelihood have been developed that can detennine these parameters from a simplified fonn 
of equation (1); 
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where Gf is the peak stress in the specimen, gqv is the specimen characteristic strength, 
and V e is known as the effective volume. 

To predict component reliability for multiaxial stress states the Batdorf theory [ 10 , 1 1 ] is 
used. Batdorf theory combines the weakest link theory with linear elastic fracture mechanics. 
It includes the calculation of the combined probability of the critical flaw being within a 
certain size range and being located and oriented so that it may cause fracture. 

Slow crack growth refers to the stable extension of a crack over time. It results from 
the combination of stress at the crack tip and chemical attack such that chemical bonds 
break and the crack tip extends. At high temperatures it can also occur from the action of 
stress and the loosening of viscous phases such that material at the crack tip displaces. 
The crack length, a, as a function of time, t, can be expressed as a power law [ 1 2 ] with 
the following form 


^P=A(x,y,z,t) K3‘"’"CF,t) (3) 

where Ki eq is the equivalent mode I stress intensity factor from the applied effective 
stress, ( x F,t) denotes a term that is a function of 'F (indicating a location (x,y,z) and crack 
orientation described by angles a, and (3 ) and t, and A and N are material parameters that 
depend on the temperature and enviromnent and hence become a function of time and 
location but not orientation. 

To take into account the time dependence of loading and material response, we 
discretize the stress history for each finite element i into short time steps (At,) during 
which the stress and material parameters are assumed to remain constant over that time 
step interval. For a specific time step j, the applied equivalent stress in element i is given 
by crieqj, the temperature Tj, the scale parameter aovj, the Weibull modulus m V j, the 
fatigue constant Byj, and the fatigue exponent Nyj. At time tk (for k time steps) the 
component survival probability, using the Batdorf approach, can be expressed as 
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B Vj = 


Aj Y 2 K IC>j Nyr2 (N V j- 2) 


(6) 


2 

The fatigue parameter B has units stress x time. Y is the crack geometry correction 
factor, Kic is the critical mode I stress intensity factor, ai eq ,k,Tmax represents the effective 
stress examined over k time steps which yields the highest fast-fracture failure probability 
[5-5], and dQ = sina da dp. Most of the terms in these equations are functions of 'F and 
t. Equation (4) was developed based on an approach of using flaw strength and 
maintaining compatibility of failure probability between discrete time steps. This 
methodology allows for the introduction of a variable Weibull modulus as a function of 
time or temperature. The effective stress Gi eq j represents an equivalent normal stress on 
the crack face caused by the combined action of the normal stress a n j and the shear stress 
Tj on the crack face, oriented normal to angles a and p - which define the orientation of 
the flaw. A mixed-mode fracture criterion and an assumed crack shape define the 

relationship for the equivalent stress. A parameter, C , known as the Shetty shear- 
sensitivity coefficient [75], is used in CARES/Z//e to weight the contribution of Tj relative 
to a n j for the equivalent stress. Shetty [75] has found this parameter to range between 
0.80 < C < 2.0. As C increases, the response becomes progressively more shear 
insensitive. The term kBvj in equation (5) is the normalized Batdorf crack density 
coefficient for volume flaws. It is used in the reliability equation for compatibility 
purposes. It insures that the multiaxial Batdorf theory collapses to the basic uniaxial 
Weibull equation (1) when a uniaxial stress state is applied for a single time step solution. 
Also, it needs to be clarified that for increased numerical accuracy CAKES/Life further 
subdivides each of the elements into their constituent Gaussian integration points and 
performs calculations in equations (4,5,6) based on these values. 

Repeated Block (Cyclic) Loading — Repeated block loading refers to a single cycle or 
a group of cycles (a block) that are applied repeatedly. Using equation (4) if k number of 
time steps describes the first load block, then for Z number of load blocks equation (4) 
requires k-Z number of time steps in order to calculate a reliability solution. Obviously, 
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this quickly becomes computationally intensive when reliability solutions for large block 
(or cycle) counts are desired. As an alternative the authors have developed a more 
computationally efficient means to perform this task [5]. This approximation method 
allows for a tradeoff between solution accuracy and numerical efficiency. In its most 
computationally efficient form this equation is expressed as 


P 5F (Z^) = exp{-2f 1 [j[...[[( 
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where it can be seen that only k computation steps are required to solve the equation 
regardless of the number of load blocks Z. 

Combined Fatigue and SCG Behavior — Equation (7) for component reliability is 
based on the power law shown in equation (3). This equation is affected by the shape of 
the cyclic wave form and is predicted to be less damaging than static loading - when the 
static load is the peak stress of the wave form applied over the same service time as the 
cyclic load. Glasses seem to follow this behavior well, while polycrystalline materials 
may show more complex behavior under cyclic loading. Damage from cyclic loading 
can be caused by a variety of effects, such as debris wedging or the degradation of 
bridging ligaments, but essentially it is based on the accumulation of some type of 
irreversible damage that tends to enhance the crack growth. To empirically account for 
cyclic effects the CARES/LIFE program has implemented the Walker law [14], which 
traditionally has been used for metal fatigue, superimposed with the power law [75] 


^pl = A lv (x,y,z,t)K kq CF,t) 


N v (x,y,z,t) 
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( 8 ) 


where Ai, A 2 , N, and Q are material constants which depend on temperature, K t ( v F,t) is 
mode I equivalent stress-intensity factor at location/orientation T\ K leq cycliCm x ( v h, t) is the 

maximum mode I equivalent stress-intensity factor (i.e., the maximum value of Ki eq over 
the cycle at location/orientation l F), AK lcq ( v F,t) is the range of mode I equivalent stress- 

intensity factor at 'F, and f c is the frequency. A methodology for estimating these 
parameters from cyclic fatigue specimen rupture data is described in [75]. The values of 
Ai, A 2 , N, and Q in equation (8) depend on location x,y,z and time t. We will assume that 
da(Hft)/dt is continuous at fractions of a cycle, that is, we assume da( l F,t)/dt could be 
computed for non-integer (real number) cycle counts. The advantage equation (8) has 
over the Walker law for ceramic materials is that when the range of the stress intensity 
factor AKi eq is zero, crack growth is not predicted to be zero in the presence of a static 
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load. Following the same derivation methodologies in [7,5] an equation similar to 
equation (7) can be constructed 
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where R is the R-ratio, shown below as a function of 'F and t 
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There are limitations associated with equation (9) and it’s implementation in 
CARES/ Life. The first is that the equation represents an averaging of properties over the 
cycle. There is no weighting given to whether the time step is at the peak or trough of the 
cycle — all time steps are weighted equally regarding their damage contribution. Thus, 
equation (9) is not a true thermomechanical fatigue methodology for brittle material 
transient reliability analysis. It is however usable when parameters are constant or vary little 
over the cycle. In CARES/ Life RfHft) is calculated by examining all time steps; therefore 
the transient FEA must be restricted to one cycle of loading when equation (9) is used. 

Proof Testing — Prior to placing a ceramic component in service, confidence that it 
will perform reliably is usually demonstrated through proof testing. Ideally, the boundary 
conditions applied to a component under proof testing simulate the conditions that the 
component would be subjected to in service. Additionally the proof test loads are 
designed to be appropriately greater in magnitude over a fixed time interval t p . After 
proof testing, the survived component is placed in service with greater confidence in its 
integrity and a predictable minimum service life, t m i n , where no failure should 
theoretically occur. 

The attenuated probability of failure, Pf a v, of a component surviving proof testing for 
time t p and subjected to an arbitrary (but known) service history over a time interval (t q - 
t p ) is 


, . , . . . Psiv(tq) 

P faV (tq) - l-P,SaV(tq)-l — — 

PSpv(tp) 


( 11 ) 


The tenn Ps p v(t p ) is the probability of survival of a component subjected to a proof test 
over a time interval denoted by t p . The tenn Psiv(tq) is the probability of survival of a 
component subjected to a proof test equivalent stress distribution ai eqp ( v P,t) over time 
interval t p and an in service equivalent stress distribution ai eq ( v P,t) over time interval (t q - 
t p ). The reliability of the survived component increases as the ratio of the proof test 
stress to the service stress increases. 

In the CARES/Z//e program the computation associated with equation (11) is not 
perfonned at the component level, but rather at 'P, a given x,y,z location and flaw 
orientation a, p. Therefore the proof test loading does not have to exactly simulate the 
service loading. For example, the proof test load and the service load can be applied in 
different directions. Further details of this methodology are described in [7] and the 
extension to transient reliability analysis is reasonably straightforward. 

Examples 

Example 1: Thermal Shocked Disk In Fast-Fracture 

In this example we examine the fast-fracture reliability response of a laser induced 
thennal shocked disk made of silicon nitride. The purpose is to see if we can predict the 
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strength response of the disks from rupture data of simple beams in uniaxial flexure. We 
will also compare the transient fast-facture reliability predictions of the disk versus the 
instantaneous fast-fracture reliability predictions of a given time step in order to verify 
the correctness and accuracy of the numerical algorithm. This example is derived from 
an international study involving laboratories from Gennany, Japan, and the United States 
[16\. Phase I of that study worked to develop and verify thermal up-shock techniques in 
which a disk specimen was centrally heated to fracture by an appropriate heating source, 
including laser, quartz lamp, shaped heating element, and gas torch. Phase II of the study 
involved a round-robin between the three countries testing thermal fracture stress for two 
silicon nitrides, AlliedSignaTs AS800 and Kyocera’s SN282. For this example we used 
disk results from the Siemens AG organization, as provided by Rettig [77], for the SN282 
material. These disks were tested using the laser irradiation technique described by 
Kirchhoff [75] and Rettig [79], We also used three-point flexure bar data that was 
provided by Ferber [2 O']. 

Thin disks 20 mm in diameter and 0.3 mm thick were centrally heated by a 800 W 
laser working in continuous wave mode. A large centrally heated area and a steep 
temperature gradient near the edge was created, which yielded high tensile stresses near 
the edge. The specimens were rapidly heated so that fracture would occur in less than a 
second. The temperature versus time response across the disk was measured with a fast 
scanning pyrometer. The heating time and the thinness of the disk were chosen such that 
through-thickness temperature gradients (and hence bending stresses) were negligible. 
Further details regarding the experimental set-up are found in [7 6, 1 8, 1 9]. 

A total of 15 disks were fractured. The time of rupture and the radius corresponding 
to the location of fracture was recorded for each specimen. Fracture stresses were 
computed by the study participants using the temperature profile at the instant of fracture, 
temperature dependent elastic modulus and thermal expansion coefficient, and integral 
equations from standard elasticity theory. Some of the disks were cut into three-point 
flexure specimens in order to independently quantify the strength characteristics. 
Fourteen of these flexure specimens were tested at room temperature with average size of 
0.30 mm thick, 3.25 mm width, 15.00 mm long with a support span of 9.44 mm. The 
specimen edges were not beveled. 

In order to perform the CARES/Zz/e reliability analysis an ANSYS finite element 
model of the disk was prepared. The model comprised a 90° slice of the disk and spanned 
Vi the thickness (1/8 of the disk is modeled). Solid elements were used in the model. The 
disk was not constrained (other than to prevent rigid body motion) and thus was freely 
allowed to expand. Temperature dependent thennal expansion and elastic modulus values 
from [16] were used. The thermal loading profiles versus time for two specimens 
(designated as #3 and #9) as supplied by Rettig [77] were used for the thermal loading. 
The temperatures were assumed constant through the thickness. Specimen #3 had the 
highest time to failure of 0.65 seconds and, hence, the highest reported fracture strength 
of 430 MPa [16]. Specimen #9 had an intennediate fracture strength of 340 MPa [16]. 
Figure 1 (a) shows the transient thermal profile for specimen #3. Figure 1 (b) shows the 
specimen #3 transient stress analysis results from the FEA. 

The figure shows the tangential (circumferential) stresses versus distance from the 
disk center. The tangential stresses are compressive near the disk center and become 
tensile towards the disk edge. The radial stresses are always compressive and are not 
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FIG. 1 — Disk specimen #3 - (a) Transient thermal profile, and (b) Transient tangential 
stress profile. Time steps range from 0. 0 to 0. 65 seconds. Not all time steps are shown. 


shown. The FEA analysis for specimen #3 consisted of 27 time steps ranging from 0.0 to 
0.65 seconds, while the analysis for specimen #9 consisted of 15 time steps ranging from 
0.0 to 0.35 seconds. The FEA transient tangential stress results compared very well with 
the numerical calculations performed by the Siemens group. 

For the CARES/Zz/e reliability analysis we used the Weibull parameters obtained 
from the three-point flexure bars to predict the strength response of the disks. Utilizing 
equation (2) with maximum likelihood analysis and assuming volume flaws we obtained 
a Weibull modulus m v = 1 1.96, a characteristic strength aev = 612.7 MPa, and a Weibull 
scale parameter aov = 453.8 MPa-mm 3/m for the flexure bars tested at room temperature. 
Using these values of m v and aov and the results of the FEA with equation (4), the 
transient reliability response of the disk was calculated. For fast-fracture transient 
analysis, the time increment terms At can be set to zero and only ai eq ,k,Tmax has to be 
determined. Since only one set of Weibull parameters at one temperature are used, 
c?ieq,k,Tmax in this case is simply the maximum effective stress at l F of all time steps used in 
the analysis. Material strength of ceramics (and hence Weibull and fatigue parameters) 
are known to be temperature dependent. Reference [16] shows this relationship for 
SN282, where average strength gradually lowers as temperature increases. We however, 
are not considering this dependency since only room temperature results were available 
for the bars cut from the disks. We believe this is still satisfactory since from Figure 1 (a) 
the temperatures near the edge of the disk (where fracture is most likely to occur) are 
relatively low such that deviation from room temperature values should be limited. 

Figure 2 shows the predicted failure probability versus time for specimens #3 and #9. 
These results are obtained using volume based analysis with the Batdorf multiaxial 

methodology, a Griffith crack and C =0.82. The plot shows straight line segments 
connecting the failure probability predictions for the various time steps. Each time step is 
based on analysis results from the experimentally measured temperature profile. The 
solid line shows the results of the transient analysis from equation (4), while the dotted 
line are results from fast-fracture analysis of the individual time step. Notice that the 
dotted curve occasionally shows a lower failure probability than a previous time step, 
while the solid curve for the transient analysis correctly does not show this trend. Also, 
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0.-D 0.20 0.30 0.40 0.50 0.60 0.70 

Time [seconds] 

FIG. 2 — Failure probability as a function of time for disk specimens #3 and #9. 
Solid line is transient fast-fracture prediction and dotted line is fast-fracture 
prediction for a particular time step. 



FIG. 3 — Predicted failure probability of disk versus stress, ay, using the Weibull 
parameters estimated from the 3-point flexure bar data. Solid line is for specimen #3 
and dotted line is specimen #9. Experimental rupture data also shown. 

there is a close correlation between the transient fast-fracture results and the single-time- 
step fast-fracture results, increasing confidence of the validity of the new algorithm. The 
solid curve and the dotted curve of disk #9 are virtually coincident. Disk #3 and disk #9 
truncate at different failure probabilities because each failed at a different maximum 
stress df. Another interesting observation is that disk #3 and disk #9 appear to have 
somewhat different failure probability responses versus time. For potential follow-up 
work it would be interesting to compare the other disks to determine the range of 
responses. 

Figure 3 shows the predicted failure probability response of the disk versus the 
maximum stress ar in the disk. The experimentally obtained fracture stresses are overlaid 
for comparison. The curve for the three-point flexure bar results represents a line of best 
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fit to the data as previously described (m v = 1 1.96, and aev = 612.7 MPa) that was used 
to obtain the Weibull parameters used for the disk reliability analysis (m v = 11.96, and 
ctov = 453.8 MPa-mnr m ). The solid curve for the disk represents predictions based on the 
analysis of disk #3, while the more difficult to see dotted curve is the prediction from 
disk #9. Notice that the disk #9 results truncate around Pf = 0.45 consistent with 
Figure 2. The curves for disk #3 and #9 follow nearly the same path in Figure 3, unlike 
the results shown in Figure 2. A striking observation about figure 3 is the difference in 
median strength between the disk and the three-point flexure bar. This primarily 
represents the “Weibull size-effect” — meaning that a component with a lot of volume 
under high stress will have a lower average strength than a component with a small 
amount of volume under high stress, and is a direct consequence of equation (2). 

Another interesting item worth commenting on is the significant difference between 
the Weibull modulus m v for the three-point flexure bar rupture data (m v = 1 1.96, aev = 
612.7 MPa) and that of the thermal shocked disk (m v = 6.91, aev = 345.9 MPa) 
experimental rupture data (shown in Figure 3) as detennined by CARES/Zz/e maximum 
likelihood parameter estimation. Least squares regression (using an Excel spreadsheet) 
on the CARES/Zz'/e disk predictions curves shown in Figure 3 yields a Weibull modulus 
m v = 8.72. This is an interesting result because under usual circumstances the Weibull 
modulus for the test specimen and the designed component are presumed to be the same 
given that they are sampling similar flaw populations and that the effective volume V e 
stays constant. Meaning, under usual circumstances the Weibull modulus obtained from 
regression of the predicted disk failure probability curve versus strength would be 1 1 .96, 
the same value as the three-point flexure specimen data. The fact that in this case the 
Weibull modulii between the predicted desk response curve and the flexure bar data are 
significantly different can be understood by examining Figure 1 (b). In the figure notice 
that as time increases, the amount of volume under high tensile stress significantly 
decreases. This has the effect of decreasing the effective volume V e with time, and 
through the size effect, increasing predicted failure stress af. The overall effect of this 
changing effective volume is to decrease the observed (apparent) Weibull modulus for 
the thermally shocked disks based on af. The fact that CARES/Zz/e predicts that the 
Weibull modulus obtained from the failure probability versus af curve reduces to 8.7 
compared to the flexure bar value of 12.0, compares favorably to the to the experimental 
disk result of 6.9. In other words, some of the discrepancy in the Weibull modulus 
values between the experimental disk and flexure bar results can be explained as a 
consequence of the transient thermal loads and how they influence the stress distribution 
in the disk as a function of time. The difference is not necessarily indicative of some 
error or inconsistency. The remainder of the difference can be explained as natural 
statistical variation (within 90% confidence bounds). 

The excellent correlation in Figure 3 to experimental results must be considered 
within the context of the underlying statistics, given that the effective volumes V e 
between the three-point bend bar and the disk is large and sensitive to the value chosen 
for the Weibull modulus m v . Because of the relative small number of samples tested (in 
this case 14 flexure specimens) and the large size effect, using Weibull parameters based 
on 90% confidence intervals from the three-point flexure bar data could shift the disk 
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predictions significantly to the left and right of the experimental data. In general a good 
design practice would be to avoid large size-effect scalings between specimens and 
components unless experimental data exists for both specimen and component such that 
data pooling practices can be taken advantage of to obtain a set of best-fit Weibull 
parameters. 

Example 2: Diesel Exhaust Valve In Cyclic Loading 

This example, involving a heavy-duty diesel ceramic exhaust valve [21], was selected 
to contrast failure probability predictions for power law, Walker law and a proof test 
condition. The valves were made of NT-551 silicon nitride material. Table 1 contains a 
summary of the Weibull and SCG parameters obtained from four-point flexure bars at the 
three different temperatures [22]. CARES /Life was used to obtain these parameters using 
the raw data listed in reference [22]. Data to obtain values for the Walker law parameters 
did not exist, so assumed values for Q v and A 2 /A 1 were used. 

TABLE 1 — NT551 fast fracture and SCG material properties. 


T(°C) 

m v 

oov 

(MPa-mrn 3/m ) 

ere 

(Mpa) 

N v 

By 

(MPa 2 -sec) 

Qv 

A 2 /A 1 

20 

9.4 

1054 

806 

31.6 

5.44e5 

3.2 

0.65 

700 

9.6 

773 

593 

86.5 

1.12e4 

3.2 

0.65 

850 

8.4 

790 

577 

18.5 

1.13e6 

3.2 

0.65 


Fifteen valves were engine tested without failure. These valves consisted of seven 
longitudinally machined valves and eight transversely machined valves. The transversely 
machined valves had been engine tested for 1000 hours while the longitudinally 
machined valves had been engine tested for 166 hours. These valves were subsequently 
tested in fast-fracture in order to examine their retained strength. The mode of failure for 
both valve-machining orientations was found to be volume induced. 

Since all engine tested valves failed due to volume flaws, the valve’s transient 
reliability was based on volume analysis. Figure 4 (a) shows the pressure variation as a 
function of time during a typical combustion cycle of 0.0315 seconds. The pressure is 
applied to the valve’s face and other exposed surfaces within the cylinder. The maximum 
attained pressure during the combustion cycle was estimated to be 15.85 MPa [21]. A 
445 N (100 lb) force due to spring preload is applied to the valve stem when it is in the 
open position. At the moment the valve closes an impact force of 1335 N (300 lb) is 
applied to the valve stem. In addition, thermal stresses due to the temperature distribution 
in the valve are superposed to the mechanical stresses. 

Figure 4 (b) shows the approximate mean thermal profile in the valve. Steady-state 
thermal analysis using ANSYS FEA code was conducted to compute these temperatures. 
This figure shows that the temperature is maximum near the valve face and decays 
towards the valve seat and stem. 

Transient reliability analysis (using equations (7) and (9)) was conducted by dividing 
the load history into 29 time steps. During each step, the load was assumed constant. The 
loads corresponding to these time steps were modeled within the ANSYS FEA program, 
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FIG. 4 — Loading profile: (a) Pressure variation applied to the face of a ceramic valve 
during a typical engine combustion cycle, and (b) mean thermal profile. 


which yielded the stress results for these 29 time steps (stress history). Figure 5 (a) 
highlights the first principal thermomechanical stress distribution in the valve at the 
moment of maximum applied pressure (at time step 6). From the figure it is apparent that 
the maximum stress location is at the valve radius, which is in agreement with the FEA 
results of Corum et al. [21], 
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FIG. 5 — (a) First principal stress distribution in the valve at the moment of maximum 
applied pressure (MPa), (b) failure probability versus cycles for proof test (dotted lines) and 
non-proof-test (solid lines) conditions for a static load, power law, and combined law. 
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The valve’s stress history and other relevant terms (temperature, volume, material 
properties, element number, etc.) were subsequently read into CARES/ Life. Figure 5 (b) 
shows the predicted failure probability as a function of cycles for various senarios. None 
of the tested valves failed, 8 of which were tested to 1000 hours (which corresponds to 
approximately 1.1 x 10 cycles). Thus the failure rate of the tested valves was less than 1 
out of 8 at 1000 hours of operation. From Figure 5 (b) it can be seen that in fast- fracture 
(at one cycle with no SCG) about 5 out of 100,000 valves are predicted to fail from the 
loading. With SCG and the power law about 20 out of 100,000 valves fail after 1000 
hours operation. If we conservatively assume that a worst-case load (time step 6) is 
applied over the whole engine cycle (static loading), then about 40 out of 100,000 valves 
fail after 1000 hours. Using the combined Walker — power law with the hypothetical 
parameters from Table 1 this rate increases to 70 out of 100,000 valves at 1000 hours. 
All of these rates were well below what was experimentally observed. These three 
scenarios were presented in order to contrast the predictions. The power law with cyclic 
loading predicts the least damage. The static loading scenario is predicted to be more 
damaging at double the rate of the cyclic loading. The combined Walker-power law 
predicts the most damage due to the enhanced cyclic fatigue effect. Since the effect of 
cyclic loads is unknown, a more conservative assumption is to assume static loading, 
however even this assumption can be un-conservative. 

To reduce the predicted failure rate even further, proof testing can be applied to 
prevent the weakest components (those with the highest likelihood of failing) from being 
placed into service. Figure 5 (b) also shows the predicted results for the three various 
scenarios for an applied proof test of 10,000 cycles at a factor of 1.1 of the service 
loading (dotted lines in the figure). Increasing the proof test load reduces the attenuated 
probability of failure even further. 

Example 3: Alumina In Static Fatigue — Material Properties Changing with Time 

In this example we examine how the reliability response of a vitreous bonded alumina is 
hypothetically affected by Weibull and fatigue parameters that change over time. The two 
scenarios that are presented are strictly based on curve fitting the data. We do not have the 
sufficient infonnation necessary to physically specify which parameters are really 
responsible for the change in material response with time. However, the point of this 
analysis is to present an analytical model which potentially can explain the nonlinear and 
changing rupture behavior of the material given that we know what physical parameters are 
responsible for the changing material behavior. This capability is potentially useful in 
modeling materials with changing composition (or changing physics of crack growth), 
including oxidation and crack blunting/healing phenomena. 

The data for the example was obtained from Quinn [23] and consists of rupture lives of 
alumina four-point flexure bars in static fatigue (loading at a constant stress level over time). 
The specimens had average dimensions of 2.2 mm height, 2.8 mm width, a load span of 
19.0 mm, and a support span of 38.0 mm. This data set was chosen from the literature 
because of its “non-ideal” behavior and the careful experimental technique reported by the 
author, which reduces the likelihood that equipment and measurement errors significantly 
affected the results. The rupture data for the individual specimens is shown in Figure 6. 
Testing was perfonned at 1000 °C. The report stated that very little creep deformation was 
detected, and although fractography was attempted, it did not reveal the source of the 
strength limiting flaws. 
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Two trends make this data non- ideally behaved: (1) dramatically increased scatter at 
lower applied stresses, and (2) pronounced nonlinear behavior in the stress rupture data as 
the applied stresses became lower (the strength data shifted to the right which signified that 
the lives became longer than would have been predicted if the standard power law was 
applicable). An ideally behaved data set would show a straight line trend indicating that the 
fatigue exponent is constant, also the scatter band would appear with a constant width versus 
the applied stress indicating the Weibull modulus is constant. 

A confirmatory piece of evidence that material properties were changing with time was 
obtained by performing static rupture experiments on specimens that were annealed (at no 
load) for either 1 or 24 hours at a temperature of 1050 °C. The outcome of these treatments 
(not shown herein) was that life dramatically increased versus exposure time for an applied 
stress level [23]. Quinn [23] indicated that the observed behavior is probably due to the 
partial devitrification of the material’s glassy phase, resulting in an increased viscosity of the 
phase and a material with improved creep and stress rupture behavior. Quinn cites reference 
[24] to support this, although he concedes that crack blunting, healing, or residual stress 
changes could also be operative. 

The data of Quinn defies the conventional modeling approach since the material 
properties could be changing over time. To account for this we show what happens when 
Weibull and fatigue parameters change with time for two cases; (1) changing the fatigue 
exponent N and the Weibull modulus over time, and (2) changing the Weibull parameters 
(m and ao) over time. Tables 2 and 3 show the sets of parameters chosen to demonstrate 
these scenarios. Given the fact that no established parameter estimation techniques currently 
exist for this type of non-ideally behaved data, these parameters were iteratively selected to 
fit the data. Table 2 contains a set of parameters versus time where the fatigue exponent, N, 



FIG. 6 — Affect on failure probability of changing fatigue exponent, N, with the log of time 
from Table 2 on static fatigue rupture data of an alumina at 1000 °C. 
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TABLE 2 — Weibull and fatigue parameters associated with Figure 6. 


t (seconds) 

m s 

(<^0S) 

(MPa-mm 2/m ) 

N s 

B s 

(MPa 2 -sec) 

1.6 

29.4 

156.8 

6.7 

2711.1 

31.6 

15.8 

152.7 

13.2 

9707.7 

1.0e5 

13.1 

127.3 

36.4 

2276.2 


TABLE 3 — Weibull and fatigue parameters associated with Figure 7. 


t (seconds) 

m s 

(<^0S) 

(MPa-mm 2/m ) 

N s 

B s 

(MPa 2 -sec) 

1.6 

29.4 

165.8 

6.7 

2711.1 

31.6 

7.4 

263.3 

8.0 

2395.9 

316.2 

4.5 

870.1 

9.0 

10,389.0 
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FIG. 7 — Affect on failure probability of changing Weibull modulus, m, with the log of time 
from Table 3 on static fatigue rupture data of an alumina at 1000 °C. 



NASA/TM— 2003-212498 


18 


and the Weibull modulus, m, are varied to yield an improved fit to the data. Table 3 shows a 
set of parameters where the Weibull modulus, m, is varied while the fatigue exponent, N, is 
kept relatively constant. Note that in both cases the scale parameter ao is dependent on the 
Weibull modulus, and B is dependent on the fatigue exponent as well as the Weibull 
parameters. When applying these parameter tables within a reliability analysis, the Weibull 
and fatigue parameters are linearly interpolated with the log of time within the time spans 
listed in the tables and are held constant outside of the time span. 

Because of the simplicity of the four-point specimen loading and geometry, results from 
FEA were not needed for the reliability analysis. Instead, a closed form expression for the 
effective area, A e , was used to evaluate the integral of equation (2) for a surface flaw failure 
mode [25]. An effective area of 58.0 mm 2 was calculated based on a Weibull modulus 
value of 7.7. This value was chosen because sensitivity analysis using equation (2) 
indicated that the maximum error in strength for a given failure probability would be 1% or 
less for Weibull modulus values ranging between 5.0 and 30.0. This level of error is 
negligible for the illustrative purposes of this example. To model the transient material 
response 10 time steps are used per decade of the log of the time. Hence 10 time steps are 
used between 1 and 10 seconds, while 20 time steps are used between 1 and 100 seconds. 
The time steps are log increments of time — that is, they would appear as equally spaced 
increments in Figure 6 with 10 steps per decade. 

Figure 6 shows the predicted 10%, 50%, and 90% failure probability iso-lines for the 
parameters shown in Table 2, while Figure 7 shows 1%, 10%, 30%, 50%, 70%, 90%, and 
99% failure probability iso-lines for the parameters shown in Table 3. Neither plot 
represents an optimized set of fitted parameters. As previously stated no established 
parameter estimation techniques currently exist for this type of non-ideally behaved data. 
Parameters were obtained by segregating portions of the data and perfonning parameter 
estimation as described in [7] for the four-point bending bar specimen as well as further 
refinements by trial-and-error. Figure 6 shows that the curvature in the data can be captured 
by changing the fatigue exponent, N, however, accounting for the change in scatter still 
required modifying the Weibull modulus, m. Figure 7 is interesting because a satisfactory 
fit to the data can be obtained primarily by changing the Weibull modulus and scale 
parameter only (B is changing mainly in response to these parameters). The “kink” shown 
in the failure probability iso-lines between 10 and 100 seconds is not purposely modeled — 
that is, we have not selected parameters values specifically to obtain this response. The 
“outlier” rupture data (data that appears to be not part of, or not consistent with, the main 
body of data) at long times to failure at 70 MPa and the “outlier” data at short times to 
failure at 60 MPa are better accounted for in the model compared to Figure 6. That is, data 
that visually appears as “outliers” are actually consistent with the flaw population failure 
probability response — in other words they are predicted to be there. 

The modeling assumptions that produced the results in Figure 7 reasonably capture all of 
the trends in the data as well as or better than the approach used in Figure 6. Also, both 
approaches required changing the Weibull modulus to account for the range of scatter. A 
changing Weibull modulus could be indicative of new flaw generation, or it could be 
indicative of changes of the physics of crack growth associated with R-curve behavior. On 
the other hand a changing fatigue exponent with time is consistent with crack 
blunting/healing phenomena, where the crack velocity relationship in equation (5) is 
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changing with time. In this case further tests on the material would be needed to understand 
the underlying mechanism that was driving the SCG behavior. 

Conclusions 

A methodology for computing the transient reliability in ceramic components 
subjected to fluctuating thermomechanical loading was developed and incorporated into 
the CARES/Z//e code. This enables CARES/Z//e to be used to predict component 
reliability for situations such as thennal shock, startup and shutdown conditions in heat 
engines, and cyclic loading. The methodology accounts for varying material response, 
whether due to temperature or environmental changes by allowing Weibull and fatigue 
parameters to vary over the loading history. Examples demonstrating the viability of the 
technique for fast-fracture, cyclic loading, and proof testing, were presented. 
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